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The spin transfer torque in all-metal dual spin valve, in which two antiparallelly aligned pinned 
ferromagnetic layers are on the two sides of a free ferromagnetic layer with two thin nonmagnetic 
spacers in between, is studied in the ballistic regime. It is argued that, similar to the results 
in the diffusion regime, the spin transfer torque is dramatically enhanced in comparison to that 
in a conventional spin valve although no spin accumulation exists at the magnetic-nonmagnetic 
interfaces. Within the Slonczewski's approach, an analytical expression of the torque on the free 
magnetic layer is obtained, which may serve as a theoretical model for the micromagnetic simulation 
of the spin dynamics in dual spin valve. Depending on the orientation of free layer and the degree of 
electron polarization, the spin transfer torque enhancement could be tens times. The general cases 
when transmission and reflection probabilities of free layer are different from zero or one are also 
numerically calculated. 

PACS numbers: 72.25.Ba, 75.60.Jk, 85.70.Ay, 85.75.-d 



I. INTRODUCTION 

Current induced magnetization reversal of magnetic 
multilayers has attracted much research interest due to 
its rich physics and potential applications in spintronic 
devicesji"— Spin valve, consisting of two ferromagnetic 
layers and one nonmagnetic spacer in between, is one 
of such multilayer structures. In a spin valve, one of 
the magnetic layers, acting as a spin polarizer, is thick 
so that conducting electrons are polarized after pass- 
ing through it. The polarized conducting electrons will 
transfer their spin angular momentums to local magne- 
tization of the thinner free magnetic layer, resulting in 
the spin transfer torque (STT) effect first proposed by 
Slonczewski^ and Berger.^° Although the STT has many 
advantages over a magnetic field in manipulating mag- 
netization state^^ii^ a large current density is needed to 
achieve a technologically useful magnetization switching 
speedj—ii* but the associated Joule heating could affect 
device performance. Therefore, the current density re- 
duction becomes a challenging issue from the spintronics 
application viewpointr'^ 

Many efforts have been devoted to the issue, includ- 
ing using optimized time-dependent current pulse^^ pure 
spin current)^ and thermal activation. One direct 
approach is to increase the magnitude of the STT under 
a given current through unique geometry design. In 2003, 
Berger— proposed a novel magnetic multilayer architec- 
ture called dual spin valve (DSV) where the free magnetic 
layer is sandwiched between two thicker pinned magnetic 
layers with opposite magnetizations. It is predicted that 
the STT applied on the free magnetic layer should be 
much larger than that in a traditional spin valve for a 
given current in the diffusion regime.^*' The argument 
is that spins accumulate at both non- magnetic/magnetic 
interfaces of the free layer and STT is proportional to 
spin accumulations. However, it is not clear whether 
this STT enhancement can occur in DSV in the ballistic 



regime without spin accumulations. Moreover, the an- 
alytical formalism for the STT in DSV is still an open 
problem up to nowj ^^'- - These are the focuses of present 
paper. A full-quantum description of the STT, valid 
when both the mean-free path and the spin-flip relax- 
ation length are larger than the thickness of the spacers, 
is presented. Averaged (over electron phases) STT is ob- 
tained analytically within the Slonczewski's semiclassical 
approach^ when all magnetic layers are perfect polariz- 
ers. It is found that STT in DSV, depending on the 
orientation of free layer and the degree of electron polar- 
ization, can be enhanced by a factor of tens in comparison 
with that in a spin valve. The general cases of arbitrary 
transmission and reflection coefficients of free layer are 
also numerically calculated. 

This paper is organized as follows: In Sec. II, a phys- 
ical picture of STT origin in both spin valve and DSV 
is presented. Sec. Ill is the theoretical model and for- 
mulation of the electron transport through a DSV. The 
results, including the analytical expression of the STT on 
the free layer within the Slonczewski approach, are given 
in Sec. IV. Sec. V is the summary. 



II. PHYSICAL PICTURE 

It shall be useful to present first a physical picture 
about STT origin in spin valve and, in particular, its 
enhancement in DSV. Consider a spin valve, which is 
schematically shown in Fig. 1(a), consisting of a pinned 
ferromagnetic layer FM^ and a free ferromagnetic layer 
FMb separated by a nonmagnetic spacer NM. The mag- 
netizations of FM^ and FM^ are represented by unit vec- 
tors, rxiA and ms, and their saturated magnetizations. 
For the simplicity, we assume that both ferromagnetic 
layers are perfect spin filters, such that spins aligned par- 
allelly with the layer magnetization can completely trans- 
mit through the layer, while antiparallel spins arc totally 
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in Fig. 1(b). On the top of a usual spin valve (Fig. 
1(a)), an additional pinned ferromagnetic layer FMc 
with an antiparallelly aligned magnetization to FMyi, i.e., 
mc = rnA, is added so that the free layer FMb is 
now sandwiched between FM^ and FMp, separated by 
two nonmagnetic spacers NM^i and NM^, respectively. 
Similar to the case of a spin valve, right-going electrons 
transmitting through FM^i will exert a STT on FM^ 
to align rxiB with m^, as shown by the blue arrow in 
Fig. 1(b). After the electrons transmit through FM^, 
most of them will be reflected by FMc and then impinge 
again on FM^, as shown by red-dashed symbols in Fig. 
1(b). Thus, they wih exert another STT on FMb along 
— mc = m^i direction, resulting in the STT enhance- 
ment. Multiple reflections in region 2 (Fig. 1(c)) tend 
to reduce the torque but will not cancel it totally. In 
the following sections we will verify this physical picture 
through a full-quantum mechanics calculation. 



III. MODEL AND FORMULATION 



FIG. 1: (Color online) (a) Illustration of STT generated in 
spin valve. The circles with arrows illustrate the electron 
flow with spin polarization direction, (b) Schematic explana- 
tion of STT enhancement in DSV. Vb = -Efo — Ef6 is the 
applied bias. FMa,b,c are ferromagnetic layers and NMa,s 
are nonmagnetic layers, (c) Wave functions at the interfaces 
of regions (0 — 6) in DSV and the coordinates orientation. 



reflected. In this ideal closed analytical solution 

of the STT can be obtained, which will be shown later. 

Consider electrons flowing from the left to the right 
in Fig. 1(a). The right-going electrons are initially po- 
larized along direction after passing through FM^. 
They will remain their spin polarization when they im- 
pinge on FMb, as long as the spacer thickness is much 
shorter than the spin-diffusion length, which is usually 
the case in nanoscale spin valves. Because the polariza- 
tion direction of FMb is noncollinear with that of FM^i, 
an electron polarized along m.A is the superposition of 
two states along niB and — niB, so that the component 
along niB can transmit through FMb while that along 
— niB is totally reflected since FMb is a perfect polarizer. 
Thus, there will be a net angular momentum transfer, 
perpendicular to niB, from impinged electrons to FMb, 
resulting in a torque on FMb to align its magnetization 
toward niyi, as shown by the blue arrow in Fig. 1(a). 
This is the origin of STT in spin valve. It should be 
pointed out that the subsequent multiple reflections of 
electrons within the NM spacer by the two FM/NM in- 
terfaces will reduce the STT, since the reflected electrons 
(dashed symbols in Fig. 1(a)) from FMb to FM^ and 
back to FMb will exert a torque along —niA- However, 
the net torque will not be zero since the reflected elec- 
trons have smaller flux than that of the original injected 
electrons. 

Let us now examine the spin transfer in a DSV shown 



Single charge and spin transport theory in magnetic 
multilayers was well developed and approaches vary from 
classical Valet-Fert theory in the diffusive regime ma- 
trix Boltzmann equation formalism, -^ to full-quantum 
mechanical treatments. -i^"— In the present paper, we 
will adopt a full-quantum mechanical description called 
scattering matrix method^^"— that is valid for ballis- 
tic transmission. We assume that interfaces are flat 
and clean and all spin-flip processes are negligible so 
that momentum k is a good quantum number in each 
layer. Wavef unction at the interface of layers labeled by 
n = 1 — 6 as shown in Fig. 1(c) can be written as a 
two-component spinor multiplied by spatial plane wave, 

V'«,o(^)=f''"^"'^V"^'^^-°' (1) 

where a denotes propagation directions, a = L for left- 
ward and a = R for rightward. is the a;— component 
of wave vector. Xn,a,^U) spin-up (-down) probabil- 

ity amplitude, z— axis is along uia- niB is specified by a 
polar angle 6 and an azimuthal angle (j). 

Incoming and outgoing spinor states in each region 
are connected to each other by a scattering matrix<^"— 
In region n (= 1,2,3,4,5), ipn,L and ipn+i.R. are outgo- 
ing spinors while ipn,R and ipn+i,L are incoming spinors. 
They relate to each other by a scattering matrix 5„ , 

Note that i^a^L — since electrons flow from the left 
reservoir to the right one. is a 4 x 4 matrix and can 
be expressed by transmission and reflection coefficients 
of each scattering region, 

Sn={f'A, (3) 
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where t„ (f„) with n= 1, 2, 3, 4, 5 are 2x2 transmission 
(reflection) matrices for layers FM^, NM^, FM^, NMs, 
and FMc, respectively. 

We will first treat the pinned ferromagnetic layers FMa 
and FMc as perfect polarizers. Hence ti, fi, t^, and fs 
take the following forms 



tl=f5 = 



1 




fi=h = 




1 



(4) 



Since there is no scatterings in nonmagnetic layers NM^ 
and NMs and the propagation of electrons in these layers 
accumulate dynamical phases, we have 

^2 = exp {iipA) i, ii = exp (iy-s) 7, r2 = f4 = 0, (5) 

where (pA and ips are the corresponding phase shifts and 
/ is the 2x2 unit matrix. The scattering matrix for FMb 
can be expressed by the angles 9 and </> as 



h 



R( 



)rR-'{0,^), (6) 
where R {0, (/>) is the rotation that brings z into ms, 



R = e- 



■ 2 COS : 

i± . a 

6 2 sm ^ 



-e 2 sm : 
e 2 cos ^ 



(7) 



t and r are the transmission and reflection matrices when 
niB is chosen as the quantization axis 



tu 

td 



Tu 
Td 



(8) 



where f y. t^, r,u, and are transmission and reflection 
parameters. Subscripts u and d stand for spin-up (ma- 
jority) and spin-down (minority), respectively. These pa- 
rameters are complex numbers in general. f„ = 1, = 
and r„ = 0, rd = 1 if FMs is a perfect polarizer. 

To find the STT on free layer FM_b, we need to ob- 
tain the scattering states at the two interfaces of FMb. 
The spin-dependent scattering wave functions can be ex- 
pressed in terms of incoming wave tpi^R, such as Vn,a = 
Pn,a4>i,R {n = 3, 4, 6) , where matrices Pn,a are given by 



P3,R = 







PS,L = 


hPa, 


P4,R = 




P6,R = 




with 




Q = 


[('- 



(9) 



(10) 



After some algebras, we have 

Q 

P3,R 



y/ZiZ2 
D{zi,Z2) 

/iT 



P3,L ■ 

P4,R 
Pi,L ■ 



D{zi,Z2) 
D{ZI,Z2) 



D{zuZ2) 
Z2^fz{ 



0^^(^1,Z2) 
Q^H^UZ2) 

' PVA^UZ2) 
Pi\{zx,Z2) 

Pl\{ZuZ2) 
PI\{ZUZ2) 

Q^^(^l,^2) 0" 

Q^"^ (21,^2) 0_ 

Q^^(^l,22) 0" 





D{zx,Z2) 

where z\ = exp (iflipA) , Z2 = exp (VIlpb) , and 



Q^^ {Z1,Z2) = tu cos^ 2'^^'^ ^^^^ 2 



-zi 1 rdtu cos^ -^+rutd sin^ - 



(11a) 
(lib) 

(He) 

(lid) 
(He) 



Q^HzuZ2) = =:e''<'sme 



1 



iu td 



+ {zi + Z2) {rutd - rdtu) 

+Z1Z2 (turl + tdtl - tJd - tdrl^ 



D {zi, Z2) = I - zi sin^ 2'^^'^ 2 

f 2^ . 2 0\ 

-Z2 I r„ cos -+rdsm. -\ + 

TuVd + {{Tu - Tdf - {tu - tdf^ 



(12a) 



(12b) 



Z1Z2 

-^3,^(^1' ^2) = D{zi,Z2), 
1 

2' 



^3,fl(^i>^2) = ^^le'^sinf 



Tu - Td 



+Z2 {-rl + Turd + tl- tutd) cos^ - 



+^2 {rl - TuTd -td + tutd) sin^ - 



P3 L\^i'^V = Z2 [tuCOS -+tdsm - 



(12c) 



(12d) 



(12e) 



-Z2 I r„cos^ - + rd sin^ - 



2 <j ■ 2 
+ru cos -j^+Td sm - 



(12f) 



+2:1-22 



t 2 j.2\ ■ 2 ^ , (2 j2\ 2^ 

ru[rd-td)sm - (r„ - i„) cos - 



PsJl (^1) ^2) 



^'3,^(^1,^2) 



Zl 



(12g) 



The notation X^"^ {X = Q, Ps^r, P3,l) refers to the tran- 
sition amplitude from spin-up to spin-down states. 
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IV. RESULTS AND DISCUSSIONS 

A. Charge current 

An applied bias VJ, shown in Fig. 1(b) generates a 
charge current Je and a spatially dependent spin current 
Js through the device. At zero temperature, the charge 
current reads 



Je 



je (q) 



with charge current density 

hk 



Je (q) 



m ' 



(13) 



(14) 



where q is the transverse wave vector with energy E, 
k'^ + — 2mE/H^, and e, m, H are the electron charge, 
electron mass and the Planck constant. The charge cur- 
rent density can also be written as 



e — Te (zi,Z2) 

m 



with transmission coefficient 



Te {zi,Z2) 



(15) 



(16) 



In the case that electrons propagate ballistically 
through NMa and NM^, the phase shifts in normal met- 
als are given by ipA = kxlA and ips = kxls where Ia and 
Ib are the widths of NM^ and NMs, respectively. For 
sufficiently thick (much bigger than electron Fermi wave 
length but still smaller than the spin diffusion length) 
NM layers, ipA and tps vary rapidly from state to state. 
Thus, when one sums up contributions from different 
electronic states (different kx), it is justifiable to assume 
ipA and ifiB to be random^^S and zi ~ exp (i2t^^) and 
Z2 = exp{i2ipB) are equally distributed on the unit cir- 
cle of the complex plane. ^^'^^ However, one should note 
that zi and Z2 are not independent under the ballistic as- 
sumption since Z2 — (zi)^ with p = Ib/Ia- The average 
transmission coefficient is then 



{Te) 



1 



2m 



Te {(pA,P^A)d(p/ 



TeiziA^lT) 



Z\ 



dzi, 



(17) 



where C is contour \zi\ = 1. The contour integral 
for p — 1, corresponding to the symmetric DSV 
configuration, ^^■■^^ is 



1=1 



Te (Zl) 



Zl 



; zi,i 



(18) 



where zi^i is the /— th pole of function Te (zi) /zi and s is 
the total number of poles inside the unit circle. In case 



when FMb is perfect, i.e., i„ = 1 and td = {ru = 0, 
and rd — 1), function Te (zi) /zi has only a second-order 
pole Zl = inside the unit circle. Thus, we can get the 
average transmission 



(Te) itu = 1, id = 0) = 



sm 



(19) 
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e/n 
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FIG. 2: 

angle 6 for (a) t 
ltdl' =0.17; (c) 



(Color online) Average transmission coefficient versus 
= 1 and td = 0; (b) \tuf = 0.84 and 
^ = 0.79 and \tdf = 0.49; (d) \tuf = 0.66 



tu\ 



and ltd! = 0.44; and (e) \tu\ = 0.73 and \td\ = 0.54. 

Figure 2 is the average transmission coefficient versus 
angle 9 for (a) perfect polarizers; (b) (001) interface of 
Au/Fe; (c) (001) interface of Cu/Co; (d) (110) interface 
of Cu/Co; (e) (111) interface of Cu/Co. The model pa- 
rameters for (b-e) are obtained from Ref. |32| . where they 
were extracted from the first-principles calculations. All 
transmission and reflection amplitudes are assumed to 
be real. We find that the average total transmission ap- 
proach zero as 9 goes to or tt even if td is finite for 
the minority electrons, which is different from the result 
(Fig. 3 in Ref. [28| ) in traditional spin valve. The rea- 
sons arc as follows. For 9 = 0, all electrons are polarized 
along after passing through FM^. They will be to- 
tally reflected by FMq because of mc = — niA, leading 
to zero electric current. On the other hand, for 9 = tt, all 
electrons that transmit FM^ will be completely reflected 
by FMb if it is a perfect polarizer because electron spins 
are antiparallel to ms, while electrons passing through 
FMb will maintain polarization along uia in the case of 
td ^ and they will be totally reflected by FMc because 
of mc = — m^. Hence, the average transmission in DSV 
vanishes at 6* = or tt. 

The total charge current flowing through the DSV can 
be obtained by summing Eq. psp over the transverse 
momentum q. To find an analytical expression, we will 
adopt the semiclassical Slonczewski approach.— Within 
the Stoner description of magnetism and let A be the ex- 
change energy of two spin bands of FM b , one can define 
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two Fermi wave vectors and K- for majorities and 
minorities, 



K+ = y/2mE/fi^, K_ = ^2m{E - t\) jW-. (20) 

For a ferromagnetic metal, we assume that the Fermi 
energy hes above the exchange potential, and electrons 
in NMs are ideally matched with the majority electrons 
in FM, i.e., k = Kj^. The possible momentum states 
that contribute to the current can be divided into three 
ranges.- '° 

Range a: < q < K_. Electrons of both spins in these 
states contribute to charge current Ja, 



Ja 



2e- 



h 



{2-k) m Ja 



K- 



Kl 



q^qdq 



h 



3 (27r) m 



K 



2\3/2' 



(21) 



Range h: if_ < q < Only majority spin electrons 
contribute to current J;,, 



Jb 



1 

= — e 



.0 ^- I .,Kl-q^qdq 
^ {Kl-Klf'\ 



3 (27r) m 



(22) 



Range c: K+ < q. All electrons are totally reflected, 
and there is no charge current flow, i.e., Jc = 0. 

Using the conventional definition of spin polarization 



P = 



K+ -K^ 
K+ + K. 



(23) 



where n± are the majority /minority spin densities at 
Fermi level in the FMs, the ratio Ja/Jb can be written 
as a function of polarization P, 



Ja ^ [i + py 

Jb 4P3/2 



(24) 



Note that Jb is the maximal polarized current for par- 
allel FMs configuration. Then, to get the total charge 
current, it should be multiplied by the average transmis- 
sion coefficient. Thus, the total electron current is given 



Je = Ja + {Te} Jb- 

B. Spin current and spin transfer torque 

The spin currents on two sides of FM g are 



(25) 



q 



(26) 
(27) 



with spin current densities 



hs (q) 
Us (q) 



2m 
2m 



(V'l,fl^V'3,i?. - V'Il'7i/'3,l) , (28) 
(V'Ifi^V'4,fi - ^I^l^-^a^lJ , (29) 



where a — {ax,cFy,(Tz) are Pauli matrices. It is con- 
venient to recast a in local orthogonal coordinates of 
ms X {m.B X m^) , x rriA, and ms such that 

a^aiUiB X (ms x m^) + (72na-B x niA + crams, (30) 

with 



CTi 



0-2 = — 



0-3 



— sm f 


i cos ae 


cos Oe"^ 


't' sin 





ie~"t' \ 




J 


9 sin 




6^0 _ 


cos 6 1 



(31) 



The STT on FM b is equal to the difference of the spin 
currents on both sides of the ferromagnet, 

r = Jss — J4s, (32) 

and the STT density is 

T^hs-ks- (33) 

Thus, we have 

T = ainiB X (rtiB x m^) -I- 02111^ x -I- a^rriB; (34) 
where 



2m 



•4,L 



(4 

-V'3,iO'j'(/'3,L-V'I,fl,0-iV'4,fi) , (» = 1, 2, 3) . (35) 

First of all, we can show as = because of the 
particle current conservation and the absence of spin- 
fiipping. This can be understood as follows: We first 
rotate z to ms, then spinor state V'ri.a = R{Q,4>) '^n,a 
where ipn,a is the electronic state seen along ms. Then, 
each spin density term ipl^ a'^s'ipn.a in Eq. (|35|) becomes 

^n,aR^ {0,(t))(7zR{0,(i))^n,a = 'ipl^a'^z'ipn.a, SO that Spin 

up state (parallel to ms,) and spin down state (anti- 
parallel to ms) are decoupled without mixing. In the 
absence of spin-fiipping, both spin-up and spin-down par- 
ticle currents are conserved. Thus, the STT projected 
along local magnetization uib is zero. Here we have 
used identity W {9, </>) a^R {6, (j)) = cr-,- The physical con- 
sequence is that STT can only rotate the magnetization 
without change its magnitude. 

Therefore, the STT can be divided into an in-plane 
(Slonczewski) termS t\\ — aims x [ms x wia) and an 
out-of-plane (field-like) term23, t±_ — 021115 x m^. We 
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note that the out-of-plane torque wih vanish if t„, td, ru, 
and rd are real. This is because all the spins of both 
transmitted and reflected electrons are in the same plane 
spanned by and ms under the condition, and they 
can only vary in this plane. 

Parameters ai and 02 are important to understand the 
spin dynamics in DSV.^^'^^ In Ref. HH these parameters 
are chosen to vary continuously without geometry depen- 
dence, while they are calculated in the difFusice transport 
limit in Ref. |22| . To find oi and 02 in our model, one 
needs to compute quantities 



P'1,R^ iP't.R 

(36) 



since = T„.{zi){h^kx)/m, (i = l,2). Perform the 
same averaging procedure as we did on the charge current 
early, one finds 



_ -1.6 
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FIG. 3: (Color online) Slonczewski torque Fy versus 9 for 
polarization P = 0.6 under three different conditions: (a) 
tu = 1 and td = 0; (b) = 0.99 and \tdf = 0.1; and (c) 

\tu\'^ = 0.84 and \tdf = 0.17. The unit of torque is | J^. 



^Res 



(i = l,2), (37) 



where zi^j^ is the Z^— th pole of function T^. (z\) / zi inside 
the unit circle in the complex plane with Si the corre- 
sponding total pole number. 

Since only region h contributes to the spin current r'l^^ 
one has (a^) — | (Tir.) J^, [i = 1,2) . Thus the average 
STT on the free magnetic layer is 

r =51 {(^) -JeVa.B X (ms X rXVA) + 92 {0) - JeHls X niA, 

e e 

(38) 

with scalar functions 



9. (0) = 



JalJb + {Te) ' 



(z = l,2). 



(39) 



The STT T consists of two terms. The first one is Slon- 
czewski torque rj|=(7i {9) l^J^niB x (ms x mj^) . And the 
second one is field-like torque T±=g2 {0) ^JeXn.B x ^A- 

The general analytical forms of 91 (9) and 92 [9) are 
difficult to find because of the complicated residue cal- 
culations in Eqs. (|18p and (|37p . However, they can be 
numerically evaluated for any given material with defi- 
nite material parameters td, r„, and Td- In Fig. 3, 
we present the numerical calculations of the magnitude 
of Slonczewski torque r||= —171 (9) ^Je sm9 per unit cur- 
rent versus angle 9 at polarization coefficient P = 0.6 
under three different conditions, which shows that dif- 
ferent transmission probabilities have strong impact on 
the torque and may change the sign of the torque. Sim- 
ilar sign reversal of STT is demonstrated in usual spin 
valve.— 

Here, in order to directly compare our results in DSV 
with Slonczewski's result in conventional spin valve,— let 
us consider the case of ideal FMs, i.e., tu = 1 and 
td — [fu = 0, = 1). After some algebras, we obtain 



1 

2' 



= i (zi -I- zl - 2) and T,, (zi) = -| (zi - z*). The 
averaged values are then 

(T,i)(<« = l, td = Q) 

(r,,)(i„ = i, td = o) = 0. 

Thus we get the scalar 5— functions 

9i{e) - 



(40) 
(41) 



-3 + (1 + Py I {2P^/^) - (niA ■ uibY 



92 {0) = 0, 



(42) 
(43) 



which show that only the Slonczewski torque exists in 
the ideal DSV. The absence of any layer-thickness de- 
pendence in Eq. results from the phase average 
across sufficiently thick normal metal layers. The val- 
ues of gi [9 = 0, tt) are crucial to evaluate the threshold 
current density J* needed for magnetization reversal of 
the free layer since J* cx l/gi {9 = 0,7r)ji^ 

Figure 4 is the magnitude of STT per unit current as 
a function of angle 9 for various polarization P in ideal 
DSV. The STT is symmetric against 7r/2 due to the con- 
tributions from both fixed magnetic layers, which is dif- 
ferent from the result in conventional spin valvCi ^'^^i^^ 
For polarization coefficient P < 1 it generally vanishes 
at = and 9 — tt (shown in Fig. 4). However, if P = 1, 
the torque is singular at 9 — and tt. The divergence can 
be understood mathematically and from a physics view- 
point. Mathematically, the singularities a.t 9 — and tt 
are due to the factor oi gi = — {s'm9)~^ when P = 1. 
Physically this is the consequence of perfect spin filter. 
In our model, every electron transfers its angular momen- 
tum to local magnetization when it impinges the interface 
of magnetic layer whose magnetization is not parallel to 
its spini^ Meanwhile, the STT per unit current is defined 
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0.4 ^, 0.6 

e/n 




FIG. 4: (Color online) Slonczewski torque Fy versus angle 9 
for various polarization P in the ideal case where all FMs are 
perfect polarizers. The unit of torque is - Je- 



FIG. 5: (Color online) Spin transfer torque enhancement ratio 
gi (9) / gl (9) versus 9 for various polarization P in the ideal 
case where all FMs are perfect polarizers. 



as the spin transfer per transmitted electroni^ However, 
in the case oi 9 — ot 6 — tt, perfect spin filter does not 
allow any electron transmitting through the DSV, which 
leads to a zero electron transmission and results in the 
divergence. The above argument can also be applied to 
the STT divergence at 9 = it in traditional spin valve 
in the original paper by Slonczewskii^ Nevertheless, we 
find that the total STT T = if Jfems x (ms x m^) if 
P = I, does not have such singularities. 

Finally, we would like to compare the magnitude of 
STT in our DSV with that obtained in traditional spin 
valve. In Ref. [2^ Krstajic et al. calculated the 
5— function in spin valve 



-4+(l + P)'(3 + mA-mB)/(4P3/2)' ' 

which is the same as the Slonczewski's result in Ref. [^. 

We plot the ratio gi (9) / g'l (9) as a function of angle 9 
for various polarization P in Fig. 5. One can see that the 
STT is largely enhanced when 9 is acute, but approaches 
to that in usual spin valve when 9 > tt/2. Namely, the en- 
hancement value sensitively depends on both the orienta- 
tion of free layer and the polarization degree of electrons. 
Under small tilted angle 9 and large polarization P, the 
enhancement ratio is dramatic, which will substantially 
lower the required threshold current density for magne- 
tization switching. For instance, gi/g* = 11.9 if 6* = ^ 
and P = 0.8. While the enhancement ratio approaches 
to 1 (red dash line shown in Fig. 5) when 9 is close to tt. 
The reason is that there is no difference between a DSV 
and a usual spin valve when layers FM^ and FM^ are an- 
tiparallelly aligned, since FM^ would not reflect electrons 
coming out of FM^ in such case. The results are qualita- 
tively consistent with Fuchs et aVs experiment^ which 
shows that the reduction of threshold current switching 
FMs from parallel to antiparallel, with respect to FM^i, 
is substantial, while it is only of modest size from antipar- 



allel to parallel in DSV compared to that in conventional 
spin valve. 

From the qualitative physical picture and quantum 
mechanical calculations, we conclude that STT can be 
greatly enhanced in DSV compared to that in spin 
valve structure in the ballistic regime without spin ac- 
cumulations. The findings of the physics behind and 
the analytical formula of STT in this emerging geom- 
etry should be interesting for both theoretical^ii^ and 
experimenta l^^'^^ concerns. Micromagnetic simulation 
based on our new function (Eq. (|42p ) and Landau- 
Lifshitz-Gilbert (LLG) equation'^^ will be a direction of 
future research. The behavior of STT in asymmetric 
DSV when the widths of two NMs are different, i.e., 
I A 7^ Ib, is also an interesting issue for further investi- 
gation. 

Although the advantage of STT enhancement is un- 
ambiguously demonstrated in our results, the accuracy 
of the analytical formula still needs experimental confir- 
mation. We suggest to use a recently developed technique 
called spin-transfer-driven ferromagnetic resonance (ST- 
FMR)-^^ to measure the angular dependence of the STT 
in DSV. 



V. SUMMARY 

In conclusion, we derive the STT acting on the free 
magnetic layer in a DSV structure in the ballistic regime. 
A full-quantum mechanics description of the STT is pre- 
sented, which is valid for nanoscale DSVs where both 
the electron mean-free path and the spin-flip relaxation 
length are larger than the thickness of the spacers. Us- 
ing a quasi-one-dimensional model and within the Slon- 
czewski's approach, we obtained the analytical form of 
the STT when all magnetic layers are perfect polarizers. 
Similar to the results in the diffusive regime, the STT is 
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dramatically enhanced in comparison to that in a conven- 
tional spin valve although no spin accumulation exists at 
the magnetic-nonmagnetic interfaces. Depending on the 
orientation of free magnetic layer and the polarization 
degree of electrons, the STT can be enhanced by a factor 
of a few tens. Our analytical (7— function provides a theo- 
retical base for the micromagnetic simulation of the spin 
dynamics in DSV. The general cases when transmission 
and reflection probabilities of free layer are different from 
zero or one arc also numerically calculated, which shows 
that the sign of the torque may change under different 
transmission probabilities. These results should be use- 
ful for the switching current reduction in magnetization 



reversal. 
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